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EXTENSION OF THE FOKKER- PLANCK 


EQUATION 

by 

John C. Price 
ABSTRACT 

A simple collision term valid to order 1/ln A is derived for 
small amplitude waves in a uniform plasma. This result extends 
the validity of Fokker- Planck type equations from the domain 
0 < co « co 0 < k «k . to the domain 0 < co « Aoj 0 < k < k .. 

— p’ — a — p 7 — 

w i L~/‘Z. CT 'ly/ 

I. INTRODUCTION 

It is of general interest to obtain the effects of collisions on the behavior of 
waves in a plasma. Previous efforts have been directed along two main lines. 

Some authors ’ use the linearized version of the Fokker-Planck equation to 
study the effects of collisions on waves. This is certainly correct for the study 
of low frequency («~v colUsion ) long wavelength (A-A mean fre e P ath) waves, but 
it is not clear how far outside this domain the equation may be trusted. For 
example it is obviously suspect at frequencies of the order of the plasma frequency. 



1. D. A. Tidman, Phys. of Fluids 7, 1826, 1964. 

2. G. G. Comisar, Phys. of Fluids 6, 76, 1963. 

3. M. Rosenbluth, W. M. MacDonald, and D. L. Judd, Phys. Rev. 107 , 1, 1957. 
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Other authors have used collision models^ to study the desired effects. 
These models are constructed to make the resulting analysis relatively tractable, 
but they are somewhat arbitrary, for they are not based on the reduction of a 
more exact frequency dependent collision term. 

In this paper we shall construct a reasonably accurate collision term for 
plasma waves. This collision term should prove useful in several respects: for 
direct analysis of the resulting equation for the full range of collisional effects; 
for perturbation solutions leading to a relatively simple dispersion relation for 
waves; and as a valid basis for the construction of model collision terms. 

Earlier attempts by Frieman,® Berk,^ and Silin® have involved restrictions 
which will not apply here. Guernseys' result® presumably includes the result to 
be derived here (except that Guernsey considers waves in an equilibrium plasma) , 
but his analysis is difficult and his result so complicated that the present simpler 
(and less accurate) derivation seems worthwhile. 

In section II the equations of the BBGKY hierarchy are simplified to form 
the basis of the analysis. In section HI the collision term is calculated and the 
kinetic equation is obtained. In section IV we discuss the validity of the collision 
term for plasma waves. 


4. E. P. Gross and M. Krook, Phys. Rev. 102 , 593, 1956. 

5. R. L. Liboff, Phys. of Fluids, 5, 963, 1962. 

6. E. Frieman, Journal of Math. Physics, 4,410, 1963. 

7. H. Berk, Phys. of Fluids, JL7, 257, 1964. 

8. V. Silin, JETP, U, 1277, 1960. 

9. R. Guernsey, Phys of Fluids, 5, 322, 1962. 
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II. THE SIMPLIFIED EQUATIONS 


We begin with the usual equations of the BBGKY^ 
tribution f and the pair correlation function g. With jr 




n di f 
m ^Vj 'J 


dx 2 Vi 



for the one particle dis- 
j , Vj j = X t , etc, we have 

f(x ; , t) 


n (j 
m <9vj 


J 


dX.V, 



g(x,,x 2 ,t) 
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{it * -b -V,) 6(x,,x 2 , t) 


— V 

m 1 


l r i " r 2 


(A ■ A) 


*(*■• ‘) f ( X 2- *) + 


jdx 3 f(x 3 , t)v_ 


e(x 1 , X 2 , t) + m J*dx 


(X, . t) 


3 (9v, 


V. 


r l ~ r 3 l 


r “ r 

1 1 3 1 


• TvT + 0—2) 


g(X 2 , x 3 , t) 


+ ( 1 — 2 ) 


( 2 ) 


In equation 2 the terms involving three particle correlations have been 
neglected, as well as the terms involving g which are significant primarily when 
| r - r 2 | is small. The neglect of these latter terms forces the usual short 
range cutoff in the r 2 integration in 1. 


10. D. C. Montgomery and D. A. Tidman, Plasma Kinetic Theory , (McGraw-Hill, 
New York, 1964), p. 44. 
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The neglect of the term in g of equation 1 leads to the familiar Vlasov 
equation; it is this "small" term which we wish to calculate. To do this we must 
solve 2 to find g in terms of f, then insert the result in 1. This problem (in the 
linearized approximation, which we shall use here) has been solved by Guernsey, 
but the analysis and result are formidable. We shall obtain a crude (order 1/lnA) 
result after making further approximations on equation 2. 


We first estimate the magnitude of the various terms of 2, then neglect the 
less important ones. Since we do not commit outselves to a value for the wave 
frequency the term d/dt may be large or small— in general we must keep it. We 
estimate velocities and velocity gradients by the thermal velocity V th . 
Because Jfdv = 1, we take f = l/V t 3 in a volume V t 3 h . A more exact 
analysis shows that two particle correlations are cutoff (decrease rapidly) 
for | r x - r 2 1 greater than the Debye length \ d . Thus we estimate 
distances for g by a length L where k d > L > r 0 , with r 0 the closest approach 
distance. Of course g also varies due to the variation of f. Thus we must have 
L much less than the scale length of f, which restricts our analysis to waves 
with a wavelength of the order of or greater than \ d . 


Inserting these estimates in 2, we find 


d_ 

dt 


' th 
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* th 
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We collect and cancel factors, and use the fact that 



'V/ 




V 


2 

th 



to find 



1 qf (J_ _j\ 

m L 2 Wh • V t ’ h j ‘ 



+ ( 1 — 2 ) 
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'th 
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'_ L _\ 2 

\kd) 


( 1 — 2 ) 


g (4) 


Since L < we naay drop the right hand terms in g compared to those on the 

£? 

left. The resulting equation for g is sometimes called the weak coupling equation. 



+ V. 


V. +v. 



s v , 




A,) ‘) f(*>. *)< 5 > 


As a result of these approximations we must supply both long and short range 
cutoffs to the r 2 integration in equation 1. The resulting collision term will be 
valid to order 1/lnA. 

We now linearize f and g about a uniform stationary state: f = f 0 (v x ) + f 1 { , 
v x , t); g = g°( |r, - r 2 |, v p v 2 ) + g 1 (r x , r 2 , , v 2> t). We treat superscript 

one terms as small, and discard products of such terms. The neglect of time 
dependence for f 0 and g° implies that then actual time dependence should be 
slower than that of f 1 and g 1 . Thus for very low frequency waves f °must be 


6. E. Frieman, Journal of Math. Physics, 4, 410, 1963. 
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the equilibrium (Maxwellian) distribution. For waves below the collision fre- 
quency v c , the linearized theory is inappropriate and the nonlinear Fokker- 
Planck equation must be used. In keeping with the notation of the Vlasov equation 
we set E 1 = - nq J dX 2 V x 1/ | r x - r 2 1 f 1 (X 2 , t), and find the equations for f 1 and 
g 1 - 



f‘(X,, t) 

+ F £l ( r i- 

{jt + V 1 • V 1 

'V’J 

l g^Xj, x 2 



<9f° 

dv l 
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m d \ 1 
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dx 2 V, 
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s v , 
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('(X,, t) f»(» 2 ) + f»(v,) f‘(x,, t) 


(7) 


m. THE COLLISION TERM FOR PLASMA WAVES 
A. SOLUTION FOR g 1 IN TERMS OF f 1 


We consider wave solutions for f 1 of the form f 1 (r x , v x , t) = f 1 (k, v x , a>) e ik ’ r “ lWt 
where co is in the upper half u> plane, or else analytically continued^ into the 
lower half plane. Equation 6 may be written 


(-iw + ik • Vj) f 1 (k, Vj 
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E 1 (k, u) • 
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nq 2 d f f dK iK 1( 
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CO 


-&T ' J ( k> v i’ w ) (8) 


11. L. Landau, J. Phys. (U.S.S.R.) H), 25, 1946. 
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Appropriate Fourier analysis of 7 yields 


+ i (k - k) • Vj + iK • v 2 ^ g l (k - K, K, v x , v 2 ,o>) 

47riq 2 ( d d \ f K , , 0 K-k 

' f^k, v,, w ) f°(v 2 ) + 
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( 9 ) 


We insert the result for g 1 into 8 and find 



2nq 4 P f dKK 1 
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00 


) < 10 > 


B. EVALUATION OF THE K INTEGRAL 

12 

In order to evaluate the K integral in (10) we apply the Plemelj formula 
to the quantity [a> - k • v x + K • ( v x - v 2 )] noting that co is analytically 
continued from the upper half plane. 

'co-k • v x +K 1 (v x -v 2 ) = a) - k • v x + K • (v x '-'v 2 ) - i^[o J -k-v x+ K- (v, - v 2 )] 

( 11 ) 

For co « k 0 V th ~ Aftj p the principal part is odd for large K, and the re- 
sulting integral may be neglected compared to that over the delta function. 

12. D. C. Montgomery and D. A. Tidman, Ibid, p. 56. 
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The order l/ln A error is comparable to that caused by the arbitrary 
cutoff. 


The resulting collision term may be written 


J(k, Vj , oj) 




X 


[f‘(k, V,, a) f°(v 2 ) +f°(v 1 ) f'(k, 


CO 


)(f2) 


where Q is given by 


Q(a ; b) 


1 ^ 

77 


r I K | =k 0 = 27r/r 0 
•' | K | =k d = 2fr/Ad 


=JS(a-K- b)dK 

iK. 


(13) 


Here we have neglected k in the factor (K - k )/(K - k) 2 of equation 10. We 
could also displace the origin of the K integral over this factor by k, then drop 
k. This leads to a symmetric collision term, but the difference from (12) is 
not significant to within the accuracy of the expression. 

To evaluate Q we observe that the only vector available in (13) is b , so 
Q must have the form 


Q 


AI + Bbb 


(14) 



where I is the unit tensor. We dot (13) and 14 twice with b to produce 


Ab 2 + Bb 4 


l fdK 

^Jk 4 


(K • b ) 2 8(a -K • b) 




0 


( 15 ) 


Again the error is order 1/ln A compared to the term we keep. Using (15) 
(B = - A/b 2 ), and taking the trace of (13) and (14), we find 


A 



dk 

K 2 


S(a - K • b) 


(16) 


The integral is elementary when expressed in cylindrical coordinates with axis 
parallel to b . We find for Q 


Q(a; b) - In Al 
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(17) 


The generalization from single species to a multiple species plasma is 
straightforward. We use Greek subscripts to indicate particle species, and 
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combine 8, 12 and 17 to find the kinetic equation for species m- 


(-ia; + ik • v x ) 


f 1 (k, 
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(18) 


IV. CONCLUSION 

It is probably not necessary to emphasize the fact that the collision term of 
equation 18 is approximate. Clearly "order l/ln A" is a better expression for 
the error than is "l/ln A". On the other hand even the thorough solution of the 
problem by Guernsey suffers the defect of the short range (large K) divergence, 
so that the more complicated result is not formally better than that given here. 
The collision term of (18) has the virtue of simplicity. 

This collision term differs from the usual Fokker- Planck term'* in that the 
coulomb logarithm has frequency and wavenumber dependence. This de- 
pendence involves both the relative velocity of the particles | v x - v 2 | and the 
doppler shifted frequency \co - k • v|. We establish a rough criterion for the 
significance of this dependence by requiring that most particles | v| < V th be 
within the region where the logarithm is decreased for In A. This yields 
o.i % kd V th % o p . Due to the slow (logarithmic) variation we may state that 
the Fokker- Planck collision term is modified significantly at frequencies above 
the plasma frequency. 

3. M. Rosenbluth, W. M. MacDonald, and D. L. Judd, Phys. Rev. 107 , 1, 1957. 
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